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Abstract 



We derive exact expressions for the excess number of clusters b and the 
excess cumulants bn of a related quantity at the 2-D percolation point. High- 
accuracy computer simulations are in accord with our predictions. 6 is a finite- 
size correction to the Temperley-Lieb or Baxter- Temperley- Ashley formula for 
the number of clusters per site Uc in the infinite system limit; the 6„ correct 
bulk cumulants. b and bn are universal, and thus depend only on the system's 
shape. Higher-order corrections show no apparent dependence on fractional 
powers of the system size. 
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Percolation is perhaps the simplest non-trivial model in statistical mechanics. A broad 
array of techniques have been brought to bear on it by researchers in a variety of disciplines, 
and it has continued as an active research area for many years (for reviews, see In 
this work, we restrict ourselves to the case of two-dimensional systems at the percolation 
point Pc- Even this limited subject has been and continues to be studied via renormalization 
group ( and references therein), conformal field theory Coulomb gas methods [^H, 

computer simulation ( [P,|lOl and references therein), as an example of supersymmetry ||Tl| , 
and using rigorous mathematical methods (see [IH and references below). 

In this Letter we consider the average number of clusters (connected graphs) (Nc) on a 
lattice of N sites and in addition the cumulants C„ of a related quantity, Nc + Nb/2, where 
Nb is the number of bonds. We derive new exact results for the excess or finite-size correction 
term b to the former and 6„ to the latter. By use of conformal field theory and Coulomb gas 
methods explicit expressions for both quantities on a torus (including the cylinder limit), 
or a rectangle (with conformally invariant boundary conditions) are determined. These 
expressions depend only on the aspect ratio of the geometry in question. They are universal, 
and therefore apply to any particular realization, e.g. lattice and percolation type. For bond 
percolation on a square (triangular) lattice, the formula for Uc, the number of clusters per 
site in the infinite system limit, is known from the work of Temper ley and Lieb (Baxter, 
Temperley and Ashley [Q)- Thus b represents the finite-size correction to either of these 
results. However, being universal, its validity is more general and applies to any system 
at criticality. For the torus or cylinder, where b is the leading finite-size correction, our 
formulas agree extremely well with the results of computer simulations. 

Our work was motivated by the results of |jlO[ , who found numerically that b for critical 



percolation in large systems appears to be a universal quantity, dependent upon the aspect 
ratio r but independent of the type of lattice or percolation being considered. Their values 
of b{r) are plotted in Fig. 1 below, while Table I presents new results for the infinite aspect- 
ratio limit, where the torus becomes a cylinder and b{r) ^ br. Here we find the identical 
value b = 0.3608(1) for bond percolation on both square and triangular lattices, which 



demonstrates the universality of this quantity, and agrees with the conformal prediction to 
all figures. 

Hu and Lin have studied the universality of a related but difi'erent quantity, the number 
of percolating clusters, via Monte Carlo simulation [|T^. Some recent theoretical results for 
this quantity are given by Cardy [|T^]. 

We also report (Table I) new results for a few excess cumulants 6„ of the quantity 
Nc + Nb/'2., specifically 62, which corrects the fiuctuation, 63, which enters in the excess 
skewness and 64, which is related to the excess kurtosis. For the square lattice, the bulk 



value of the fiuctuation is given by Temperley and Lieb [|T2|. For percolation on a torus, 
we find non-monotonic behavior for 62 as a function of aspect ratio. Good agreement with 
computer results for 62, and 64 is obtained. 

New computer simulations were carried out to investigate the behavior of {Nc) and C„, 
n = 2, 3, 4, in the cylindrical limit (high aspect ratio torus). We used tori of dimensions LxL' 
with L'/L > 32 and L = 4, 5, 6, 7, 8, 10, 12, 16, with additional runs at L = 32, 64, 128, 256 
where we measured (Nc) only. About 10^^ lattice points were generated for each lattice 
size, requiring a total of several months of workstation computer time. Such extensive work 
was necessary to pinpoint the higher-order corrections to b and the higher cumulants. With 
these aspect ratio > 32, it was verified that no cluster percolated the long way around the 
torus, so infinite cylinders were effectively represented. Most of our calculations were done 
for bond percolation on a square lattice. We also measured (Nc) for bond percolation on 
the triangular lattice, which can also be wrapped around to make a simple cylinder, in order 
to check the universality hypothesis. In addition, new simulations for LxL systems were 
carried out. 

Our results solidify the identification of critical percolation as a realization of conformal 
field theory. This connection was first noticed some years ago 0] and made particularly 
explicit by Cardy [Q, who studied the crossing probabilities via correlation functions of 
boundary operators. Here we operate at a simpler level, employing the partition function. 
Thus our results may suggest a way to a deeper understanding of this identification — an 



independent derivation of b, for instance, would be very interesting. There is in fact recent 
progress of this type for the crossing probabihty in Voronoi percolation 

At the critical point, the partition function of the Q-state Potts model, for Q real, may 
be expressed in random cluster form (see 1171) 



Z= Q^^^+^^^Z' (1) 

graphs 

where the sum extends over all graphs (possible placements of Nb bonds on the edges of a 
lattice), Nc is the number of clusters (connected vertices, including isolated points), and the 
coupling has been set to unity. The transition is second-order for < Q < 4, with critical 
temperature given by e^" — 1 = Q^^^. Here we make use of the Q = 1 case, where the graphs 
are equally weighted, so that Z describes percolation |]TB[|. Similar results may be obtained 



for the Ising {Q = 2) and other universahty classes; these will be reported elsewhere. 
The cumulants C„ of Nc + Nb/2 follow immediately. 

Note that Ci = {Nc + Nb/2), C2 gives the corresponding fluctuation, C3 enters in the 
skewness, etc. 

At the critical point, Z is supposed to factorize 

Z = Z-Z^ (3) 

where Z is non-universal, depending on the lattice type, boundary conditions, etc., while Zu 
encodes the universal information. The corresponding universal term F„ = — In Z^ is the 
finite-size correction to the free energy; it also generates the excess quantities studied here. 
For a given geometry, it depends only on the central charge 
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7r-cos-i(VQ/2) 



(4) 
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and the shape. 

The simplest case is an infinite cylinder of width (circumference) L. Here, if we consider 
a length L' of the cylinder, Fu is given by |ll9| , |20[ 



For an infinite strip with edges (and conformally invariant boundary conditions), the formula 
is the same, except that the 6 is replaced by 24 [|^. 

Now for percolation (Q = 1, c = 0) on a cylinder (or torus) the graphs are equally 
weighted so the number of bonds is exactly proportional to the number of lattice sites. Thus 
there is no finite-size correction to {Nb). The above then leads to Ci — {Nb)/2 = {Nc) = 
n^LV + hL'/L {V » L), with 

6=^ = 0.360 844... (6) 
24 ^ ^ 

The agreement with computer simulation results is excellent, as shown in Table I. For the 

strip (with edges), the r.h.s. of is divided by four. 

For a torus, the universal factor (for arbitrary Q) may be evaluated via the work of 

di Francesco, Saleur and Zuber 0], who employ a Coulomb gas formulation. This results in 

an expansion in powers of the parameter q = e"^'^^ , where r = L'/L (length/width) of the 

torus. For percolation, the above procedure then leads to (Nc) = ricLL' + b, with 

b = ^r + q"/' (-1 + 2v^r) + q' (-1 + v^r) 

+ . . . + g5/48 ^ 2g53/48 _ ^23/16 ^ ^77/48 ^ _ _ _ ^ (7) 

keeping all terms through 0{q^). The terms proportional to r appear because the powers 
of q in Zu depend on Q, and thus contribute to the derivative. They include the leading 
result @ for r —>■ oo, which arises from the derivative of the partition function of the unit 
operator. Note that the term correcting Z^ for "cross-topology" clusters is proportional 
to Q — 1 (see (4.8)), and thus does not contribute this way, since Q = I here. Fig. 1 
compares (|^) with computer simulations. Excellent agreement is found again. 

By symmetry, b is the same for r and 1/r, although this is not evident from the few 
terms exhibited in (|^). A graph of b vs. r shows that the leading behavior in the cylinder 
limit (the first term on the r.h.s.) is corrected by a positive term that goes smoothly to zero 
as r increases from 1. This term is such that r = 1 is the only minimum of b. 



For a rectangle, with edges, of length L and width L', Zu is known (for arbitrary c) 



via the results of Kleban and Vassileva ||2T]]. One can therefore evaluate F„ for any Q by 
the same procedure. For percolation one has (here, the finite-size correction to {Nb) is 
independent of L' /L, and is therefore not included) 

6=^lnLL'-^ln[r?(g)r?(g)], (8) 

where t] is the Dedekind ?7-function, q = e"^'^'' with r = L' /L, and q = e~'^'^l'^ . However, we 
have not attempted computer simulations in this case (or for the strip with edges), since the 
leading finite-size correction will be due to a non-universal edge term, making the precise 
determination of h or 6„ difficult. Note that the Q-dependence of 6„, for any n, involves the 
same function of q shown in (^, by contrast to the torus. This is related to the fact that Z 
on a rectangle generally couples only to the unit operator Eq. (H) also assumes there 
is no Q-dependence of the matrix element or constant term in F„ (see for details). 



The fiuctuation in Nc + Nb/'2^ is given by C2. For percolation on a cylinder, its finite-size 
correction 62-^'/-^ follows from (|^). We obtain 

&2 = ^-t|- = 0.061513... (9) 

Computer simulations results for this quantity are illustrated along with (§) in Table I. 

The excess fiuctuation of the number of clusters Nq would also be of interest. However, 
to calculate this requires the excess for (NcNb), which is apparently not accessible. A 
similar situation occurs for the leading fiuctuations in the Temperley-Lieb theory ||10| , p!2| . 

For 63 and 64 we find, similarly from (|)-(|), 63 = 5^/3/36 - 9/87r + 27V^/6A7i^ = 
-0.043 500 .. . and 64 = 25\/3/108 - 9/ An + 27v^/167r2 - 81/6471^ = -.059 933 . . .. Good 
agreement with simulations is obtained here as well; see Table I. These values are obtained 
by linear fits to the data vs. L~^, which also yields the bulk values of C2, C3 and C4. Our 
results for C2 and C3 are in agreement with the predictions 0.039 445 . . . and 0.012 913 . . . 
that follow from []T^; we have not been able to obtain C4 either analytically or numerically 



from that theory, but find from our simulations C4 ~ 0.0028. Details will be presented 
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elsewhere. The bn results are less accurate than those for b because of the problem of 
extrapolating to infinity, and because the error grows rapidly as n increases (since higher 
moments are involved). 

The above also gives the excess fluctuations 62 for percolation on a torus. One finds 

- ^5/24 



b,= b,r + q'l'' (l - i^r) - q 

^^53/48 L _ 1^3 \ _ 4^29/24 
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53^3 9 



r + KA+... (10) 



24 47r 

keeping terms through 0{q^^^) to show the appearance of in the higher coefficients. At 
r = 1, 62 = 0.105 436 634. . . is at a local maximum (note b2{r) = b2{l/r)). As the aspect 
ratio r increases, 62 decreases to a minimum of 0.103 341 600 . . . at r 1.6177. As r increases 
further, 62 increases, with limiting (large r) behavior agreeing with (^). Such non-monotonic 
behavior with aspect ratio is very reminiscent of the finite-size correction to the specific heat 
(which is itself proportional to the excess of the fluctuation of Nb) in the Ising model critical 
region. This also exhibits non-monotonic behavior as a function of r |^ (see their Fig. 6). 
This curious feature has been explained in terms of a one-dimensional array of domain walls 
|23| ; a connection to whatever mechanism is operating at the percolation point would be 
very interesting. 

We have also considered the convergence of the excess quantities. This is equivalent to 
studying the higher order finite-size corrections to the density or cumulants. In all cases 
examined, the leading such term is proportional to L~^; however its coefficient is definitely 
not universal. Fig. 2 illustrates this for b{r = 1). For the square lattice on the cylinder, 
assuming the exact ric and b, an analysis of the higher-order corrections shows no sign of a 
fractional power of L, signifying an irrelevant singularity in the renormalization group, as 
are seen in some other problems in percolation ||^, and suggests a simple analytic series 

(Nc) b 0.180 0.69 , , 
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Eq. ([TTD fits all our measurements of (Nc) to nearly their full accuracy, about ±5 ■ 10 ^. This 
form is consistent with corrections to scaling due to the breakdown of rotational invariance 



on a square lattice 



In summary, we have shown that at the percolation point in two dimensions, the uni- 



versality of the excess cluster number b, observed to hold numerically in |]TD|], follows from 
conformal field theory. We have derived explicit expressions for b and also for the excess 
cumulants of a related quantity that are in complete agreement with simulation results. 

We acknowledge useful conversations with J. L. Cardy, F. Y. Wu, and thank C. Lorenz 
for obtaining numerical results for the triangular lattice. RZ was supported by the US 
National Science Foundation under Grant No. DMR-9520700. 
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FIGURES 

FIG. 1. Excess number of clusters b vs. aspect ratio r at the percolation point. Solid line: exact 
results; points: simulations [10]. The simulations combine several types of percolation and lattice 
types. The error is less than the size of the points. 

FIG. 2. bes{r = 1) vs. L~^. Points: simulation results for an L x L square lattice, Solid 
line: fitting curve 0.88358 + 0.18L~^ consistent with (11). 6eff is defined as {Nc) — ricL^, where 
nc = (3 • 3^/^ — 5)/2 [10,12]. The conformal prediction for 6(r = 1) from (7) is 0.883 576 .... 
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TABLES 

TABLE I. Exact and simulated results for excess cumulants at the percolation point. See text 
for analytic expressions. 



Quantity Exact Simulation 

b 0.360 844... 0.3608(1) 

62(r = l) 0.105437... 0.106 

62 0.061 513 .. . 0.068 

63 -0.043500... -0.041 

64 -0.059 933... -0.059 
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